Transportation dynamics on networks of mobile agents 



O 

(N 

Oh! 
D 
00 



Han-Xin Yang,i El Wen-Xu Wang,^ Yan-Bo Xie/ Ying-Cheng Lai,^'^ and Bing-Hong Wangi ^Q 

'Department of Modem Physics, University of Science and Technology of China, Hefei 230026, China 
^School of Electrical, Computer and Energy Engineering, 
Arizona State University, Tempe, Arizona 85287, USA 
^Department of Physics, Arizona State University, Tempe, Arizona 85287, USA 
Research Center for Complex System Science, University of Shanghai for Science 
and Technology and Shanghai Academy of System Science, Shanghai, 200093 China 
(Dated: September 6, 2011) 

Most existing works on transportation dynamics focus on networks of a fixed structure, but networks whose 
nodes are mobile have become widespread, such as cell-phone networks. We introduce a model to explore the 
basic physics of transportation on mobile networks. Of particular interest are the dependence of the throughput 
on the speed of agent movement and communication range. Our computations reveal a hierarchical dependence 
for the former while, for the latter, we find an algebraic power law between the throughput and the communica- 
tion range with an exponent determined by the speed. We develop a physical theory based on the Fokker-Planck 
equation to explain these phenomena. Our findings provide insights into complex transportation dynamics aris- 
ing commonly in natural and engineering systems. 
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I. INTRODUCTION 

Transportation processes are common in complex natural 
and engineering systems, examples of which include trans- 
mission of data packets on the Internet, public transporta- 
tion systems, migration of carbon in biosystems, and virus 
propagation in social and ecosystems. In the past decade, 
transportation dynamics have been studied extensively in the 
framework of complex networks ifH-floll , where a phenomenon 
of main interest is the transition from free flow to traffic con- 
gestion. For example, it is of both basic and practical inter- 
est to understand the effect of network structure and routing 
protocols on the emergence of congestion fuHH. Despite 
these works, relatively little attention has been paid to the role 
of individual mobility. The purpose of this paper is to ad- 
dress how this mobility affects the emergence of congestion 
in transportation dynamics. 

The issue of individual mobility has become increasingly 
fundamental due to the widespread use of ad-hoc wireless 
communication networks. The issue is also important in other 
contexts such as the emergence of cooperation among individ- 
uals ll20ll and species coexistence in cyclic competing games 
ll2Tll . Recently, some empirical data of human movements 
have been collected and analyzed Iq3l l23ll . From the stand- 
point of complex networks, when individuals (nodes, agents) 
are mobile, the edges in the network are no longer fixed, 
requiring different strategies to investigate the dynamics on 
such networks than those for networks with fixed topology. 
In this paper, we shall introduce an intuitive but physically 
reasonable model to deal with transportation dynamics on 
such mobile/non-stationary networks. In particular, we as- 
sume in our model that communication between two agents 



is possible only when their geographical distance is less than 
a pre-defined value, such as the case in wireless communica- 
tion. Information packets are transmitted from their sources 
to destinations through this scheme. To be concrete, we as- 
sume the physical region of interest is a square in the plane, 
and we focus on how the communication radius and moving 
speed may affect the transportation dynamics in terms of the 
emergence of congestion. Our main results are the following. 
Firstly, we find that congestion can occur for small commu- 
nication range, limited forwarding capability and low mobile 
velocity of agents. Secondly, the transportation throughput 
exhibits a hierarchical structure with respect to the moving 
speed and there is in fact an algebraic power law between the 
throughput and the communication radius, where the power 
exponent tends to assume a smaller value for higher moving 
speed. To explain these phenomena in a quantitative man- 
ner, we develop a physical theory based on solutions to the 
Fokker-Planck equation under initial and boundary conditions 
that are specifically suited with the transportation dynamics 
on mobile-agent networks. Besides providing insights into is- 
sues in complex dynamical systems such as contact process, 
random-walk theory, and self-organized dynamics, our results 
will have direct applications in systems of tremendous impor- 
tance such as ad-hoc communication networks i24l426ll . 

In Sec. |II] we describe our model of mobile agents in terms 
of the transportation rule and the network structure. In Sec. 
Hin we present numerical results on the order parameter, the 
critical transition point and the average hopping time. In Sec. 
IIVI a physical theory is presented to explain the numerical 
results. A brief conclusion is presented in Sec. |Vl 



II. TRANSPORTATION RULE AND NETWORK 
STRUCTURE OF MOBILE AGENTS 
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In our model, agents move on a square-shaped cell of 
size L with periodic boundary conditions. Agents change 
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FIG. 1: (Color online) Schematic illustration of transportation pro- 
cess among mobile agents. The dashed circles denote the commu- 
nication range, the arrows denote the moving directions, and each 
agent is specified by a number. Agents without data packets are in 
green (light gray) and agents holding data packets are in red (dark 
gray). In the left panel, agent 3 and 4 can communicate with each 
other, where the former holds a packet. Agent 3 delivers the packet 
to agent 4 and then all agents move, as shown in the right panel. After 
the movements, all agents randomly set new directions. There can be 
more than one agent in the communication range of any agent. The 
relevant communication networks are shown in the bottom. 




FIG. 2: (Color online) The proportion of nodes P{k) as a function 
of degree k for different communication radius a. The inset shows 
that the relative size S of the largest connected component and the 
clustering coefficient (E) as a function of communication radius a. 
The number of agents is 1200 and the size of the square region is 
L = 10. 



their directions of motion 9 as time evolves, but the moving 
speed V is the same for all agents. Initially, agents are ran- 
domly distributed on the cell. After each time step, the posi- 
tion and moving direction of an arbitrary agent i are updated 
according to 



Xi{t + 1) = Xi{t) + vcos9i{t), 



yi(t + 1) = yi{t) +vsm9i{t), 



^{t) - 



(1) 



(2) 



(3) 



where Xi{t) and yi{t) are the coordinates of the agent at time 
t, and is an A^-independent random variable uniformly dis- 
tributed in the interval [—tt. tt]. Each agent has the same com- 
munication radius a. Two agents can communicate with each 
other if the distance between them is less than a. At each time 
step, there are R packets generated in the system, with ran- 
domly chosen source and destination agents, and each agent 
can deliver at most C packet (we set C = 1 in this paper) 
toward its destination. To transport a packet, an agent per- 
forms a local search within a circle of radius a. If the packet's 
destination is found within the searched area, it will be de- 
livered directly to the destination and the packet will be re- 
moved immediately. Otherwise, the packet is forwarded to a 
randomly chosen agent in the searched area. The queue length 
of each agent is assumed to be unlimited and the first-in-first- 
out principle holds for the queue. The transportation process 
is schematically illustrated in Fig.[T] 

The communication network among the mobile agents can 
be extracted as follows. Every agent is regarded as a node 
of the network and a wireless link is established between two 



agents if their geographical distance is less than the communi- 
cation radius a. Due to the movement of agents, the network's 
structure evolves from time to time. The network evolution as 
a result of local mobility of agents is analogous to a locally 
rewiring process. As shown in Fig. [T] nodes 1 and 2 are dis- 
connected while node 3 and node 4 are connected at time t. 
At time t + 1, nodes 1 and 2 depart from each other and be- 
come disconnected while nodes 3 and 4 approach each other 
and establish a communication link. Note that the mobile pro- 
cess does not hold the same number of links at different time, 
which is different from the standard rewiring process where 
the number of links is usually fixed. 

We define an agent's degree at a specific time step as the 
number of links at that moment. Figure |2] shows that the de- 
gree distribution of networks of mobile agents exhibits the 
Poisson distribution: 



P{k) = 



-W(fc)fe 



(4) 



where k is the degree, P{k) is the proportion of nodes with de- 
gree k and (fc) is the average degree of network. As shown in 
Fig. 121 the average degree (fc) increases as the communication 
radius a increases and the peak value of P{k) decreases as a 
increases. We also investigate the relative size of the largest 
connected component and the clustering properties of the net- 
work in terms of the clustering coefficient. The relative size S 
of the largest connected component is defined as 



c _ "5*1 



(5) 



where Si and is the size of the largest connected component 
and the total network respectively. The clustering coefficient 
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FIG. 3: (Color online) (a) Order parameter r;(i?) versus R for differ- 
ent values of speed v and (b) dependence of the critical Rc on v for 
different communication radius a. The number of agents is 1200 and 
the size of the square region is L = 10. The delivery capacity C of 
each agent is set to be unity. In (a), a = 1.2 and r\{F() is obtained 
by averaging over 10^ time steps after disregarding 2 x 10* initial 
steps as transients. The results in both (a) and (b) are obtained by an 
ensemble average of 20 independent realizations. The lines in (b) are 
theoretical prediction from Eqs. l l26t and l l28b . 



FIG. 4: (Color online) (a) Critical value i?c as a function of the com- 
munication radius a for different values of v. The lines are theoretical 
predictions from Eqs. l l26b and ( 128b . (b) Average hopping time (T) 
as a function of a for different values of v. The theoretical curves are 
obtained from Eqs. J20t-(l23t. For v — 0.4 and u = 1, / = 1 is used. 
The average time is computed in the free flow state. Each data point 
is obtained by averaging over 20 independent realizations. 



Ei for node i is defined as the ratio between the number of 
edges Ci among the fcj neighbors of node i and its maximum 
possible value, ki{ki — l)/2, i.e., 



E, 



2e, 



(6) 



The average clustering coefficient (E) is the average of Ei 
over all nodes in the network. The insert of Fig. |2] shows that 
S and {E) increase as the communication radius a increases. 
In particular, when the value of a exceeds a certain value, e.g., 
0.4, high values of S and (E) is attained. We also note that the 
motion speed does not influence the statistical properties of 
the communication network. In general, the communication 
network caused by limited searching area and mobile behavior 
is of geographically local connections associated with Poisson 
distribution of node degrees and dense clustering structures. 



m. NUMERICAL RESULTS 

To characterize the throughput of a network, we exploit the 
order parameter rj introduced in Ref. fl]: 



t^oo R At 



(7) 



where ANp = Np{t + At) — Np{t), (• • •) indicates the average 
over a time window of width At, and Np{t) represents the 
total number of data packets in the whole network at time t. 
As the packet-generation rate R is increased through a critical 
value of Rc, a transition occurs from free flow to congestion. 
For R < Rc, due to the absence of congestion, there is a 
balance between the number of generated and that of removed 
packets so that (ANp) ~ 0, leading to ri{R) — 0. In contrast, 
for R > Rc, congestion occurs and packets will accumulate 
at some agents, resulting in a positive value of ri{R). The 



traffic throughput of the system can thus be characterized by 
the critical value Rc which is on average the largest number 
of generated packets that can be handled at each time without 
congestion. 

Figure [Sja) exemplifies the transition in the order param- 
eter from free flow to congestion state at some critical 
value Rc- We find that Rc depends on both the moving speed 
V and the communication radius a. Figure[3jb) shows the de- 
pendence of Rc on V for different values of a. We observe a 
hierarchical structure in the dependence. Specifically, when v 
is less or larger than some values, Rc remains unchanged at 
a low and a high value, respectively, regardless of the values 
of V. The transition between these two values of Rc is contin- 
uous. The hierarchical structure can in fact be predicted the- 
oretically in a quantitative manner (to be described). Figure 
lUa) shows the dependence of Rc on a for different values of 
V, which indicates an algebraic power law: Rc ^ a^, where (3 
is the power-law exponent. We find that the power law holds 
for a wide range of a and the exponent /3 is inversely coiTe- 
lated with v. For example, for v — 0, f3 3.5 but for large 
values of v, say t; = 5, we have /3 w 2. When a reaches 
the size of the square cell, Rc is close to as every agent al- 
ways stays in the searching range of all others and almost all 
packets can aiTive at their destinations in a single time step. 

To gain additional insights into the dependence of Rc on 
the parameters a and v so as to facilitate the development of 
a physical theory, we explore an alternative quantity, the av- 
erage hopping time (T) in the free flow state which, for a 
data packet, is defined as the number of hops from its source 
to destination. As we will see, (T) can not only be calcu- 
lated numerically, it is also amenable to theoretical analysis, 
providing key insights into the theory for Rc- Representative 
numerical results for (T) are shown in Fig.Ub). We see that 
for large v, (T) scales with a as a^^ and, as both v and a are 
increased, (T) decreases. 
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IV. THEORY 

We now present a physical theory to explain the power law 
behaviors associated with (T) and then Re- A starting point is 
to examine the limiting case of t; = 0, where (T) can be esti- 
mated analytically. In particular, assume that a particle walks 
randomly on an infinite plane. There are many holes of radius 
a on the plane. Holes form a phalanx and the distance between 
two nearby holes is L. The particle will stop walking when it 
falls in a hole. The underlying Fokker-Planck equation is 



dt 



= [AV' + Uir)]P{r,t), 



(8) 



where P(r, t) is the probability density function of a particle 
at location r at time t, A is the diffusion coefficient, U (r) is 
the potential energy, U{r) = — oo inside holes and U{r) =0 
outside holes, and P(r, t) — inside holes and P{r, t) > 
outside holes. Making use of solutions to the eigenvalue 
problem: 



[AV^ + ;7(r)]$„(r) = -A„$„(r), 



(9) 



where ^nlr) is the normalized eigenfunction and A„ is the 
corresponding eigenvalue, we can expand P{r, t) as 



P(r,0 = ^c„e-^"*$„(r), 



(10) 



where Cn = J P{r, 0)<i>„(r)(ir and the initial probability den- 
sity P(r, 0) is distributed over a region of a typical size a. The 
probability that a particle still walks at time t is: 

/oo 
P(r,i)dr = ^c„d„e-^"*, (11) 
n=l 

where dn — J $ri(r)(ir. Since the n = 1 term is dominant, 
we have 



Q{t) ^ cidie-^'\ 
which gives the average hopping time as 

1 



Since 



(T) « ^. 



(12) 



(13) 



(14) 



the infinite-plane problem can be transformed into a problem 
on torus: 

V^$i(r) = -/i^$i(r), for a <r <b, 

$i(r = a) = 0,$'i(r = fe) = 0, (15) 

where h — L/^/tt and 

$i(r) = BiJoihr) + B2No{hr). (16) 



Jq and Nq are the first-kind and the second-kind Bessel Func- 
tion, respectively, and B2 = —BiJQ{ha)/NQ{ha). The quan- 
tity h can be obtained by 



4{hb) _ %{hh) 
Jo{ha) No{ha)' 



Using 



Jo{x) 



2 4 
X x^ 

64' 



Noix) 



^[(ln|+0.5772)(l-:^ + ^) 



X 



3a^ 
128 



(17) 



(18) 



], (19) 



and combining Eqs. (|9]l and ( fTSl ), we get Ai = Ah^. For 

At = 1, ((Ar)2) = we have A = ((Ar)2)/4 = 

and Ai = Ah'^ = a^/i^/S. Finally, we obtain (T) as 



(20) 



For V > 0, P(r, t) becomes zero in the area where a hole 
moves and Q{t) decays with time under two mechanisms: dif- 
fusion at rate Ai and motion of holes at the rate A^. Thus, we 
have 



(T) 



/Ai + A„ ' 



(21) 



where / is the weighting factor (0 < / < 1) that decreases 
as V increases. Specifically, / = 1 for small v and / = for 
large v. The quantity A^ is given by 



(i) A, 



(ii) A, 



(iii) A, 



Jq r$i(r)arccos( )dr 
TT r$i {r)dr 
r$i(r)arccos( )dr 



TT r$i (r)dr 



62 



L2 



(22) 



where ~ (|vi — V2 p) = \/2v, (i) is valid for < < 
2a < 6 — a or < < b—a < 2a, (ii) is valid for 2a < Vs < 
b—a and (iii) is for Vg > b—a. For large values of v, agents are 
approximately well-mixed so that we can intuitively expect 
the average time (T) to be determined by the inverse of the 
ratio of the agent's searching area and the area of the cell: 



(T) 



(23) 



The validity of this equation is supported by the fact that the 
ratio of the two areas is equivalent to the ratio of the total 
number of agents to the number of agents within the searching 
area. The estimation of (T) for large v is consistent with the 
theoretical prediction from Eqs. ( 1211 1 and (l22]i(iii) by inserting 
/ = 0. The theoretical prediction is in good agreement with 
simulation results, as shown in Fig.Ub). 
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With the aid of Eqs. ( |20] | and (l23T l for (T), we can derive 
a power law for i?^. In a free-flow state, the number of dis- 
posed packets is the same as that of generated packets Rt in a 
time interval t. For v = 0, the number of packets rii passing 
through an agent is proportional to its degree. This yields 



R{T)th _ R{T)tki 
N{k) ' 



(24) 



where ki is the degree of i, the sum runs over all agents in the 
network, and (fc) is the average degree of the network. During 
t steps, an agent can deliver at most Cit packets. To avoid 
congestion requires rii < Cit. If all agents have the same 
delivering capacity C, the transportation dynamics is domi- 
nated by the agent with the largest number of neighbors and 
the transition point Rc can be estimated by 



Rc{T)tk„ 
N{k) 



Ct, 



(25) 



where kmax is the largest degree of the network. Thus, for 
w = 0, we have 



Rr — 



NC{k) 



(26) 



where (fc) ~ Ntto^/L'^. Since the degree distribution follows 
the Poisson distribution: P{k) = e^'^'^^ {k)^ /k\, the quantity 
kmax can thus be estimated by 



1 

N' 



(27) 



Inserting (T), (fc) and kmax into Eqs. ( |26] l. we can calculate 
Rc for low moving speed v. 

For large v, Eq. ( |26] | can also be applied but kmax — {k) 
and (T) = I? ji^a} . Hence, R^ for large v is given by 



Rc = ~jjr'^°- ■ 



(28) 



Equation ( |28] ) indicates that i?c scales with a?, which is in 
good agreement with simulation results shown in Fig.UJa). 



V. CONCLUSION 



In conclusion, we have introduced a physical model to 
study the transportation dynamics on networks of mobile 
agents, where communication among agents is confined in 
a circular area of radius a and agents move with fix speed 
V but in random directions. In general, the critical packet- 
generating rate Rc at which a transition in the transporta- 
tion dynamics from free flow to congestion occurs depends 
on both a and v, and we have provided a theory to explain 
the dependence. Our results yield physical insights into crit- 
ical technological systems such as ad-hoc wireless commu- 
nication networks. For example, the power laws for the net- 
work throughput uncovered in this paper can guide the design 
of better routing protocols for such communication networks. 
From the standpoint of basic physics, our findings are relevant 
to general dynamics in complex systems consisting of mobile 
agents, in contrast to many existing works where no such mo- 
bility is assumed. 
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